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Abstract 

We study classical dynamics in the spherical ansatz for the SU{2) gauge and Higgs fields of the 
electroweak Standard Model in the absence of fermions and the photon. With the Higgs boson 
mass equal to twice the gauge boson mass, we numerically demonstrate the existence of oscillons, 
extremely long-lived localized configurations that undergo regular oscillations in time. We have 
only seen oscillons in this reduced theory when the masses are in a two-to-one ratio. If a similar 
phenomenon were to persist in the full theory, it would suggest a preferred value for the Higgs 
mass. 
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INTRODUCTION 



In a wide range of nonlinear field theories, the existence of nontrivial static solutions has 
been well established 0,0). However, nonlinear field theories can also support "oscillons," 
also called "breathers," which are localized configurations that oscillate without dissipation. 
Exact solutions of this kind include sine-Gordon breathers and Q-balls 0. A number of 
numerical and approximate analyses have established the existence of long-lived breathers 
in theories with even less structure, such as 4 theory in one dimension [5||, 4 theory in 
three 0,0,0 and higher 0] dimensions, higher dimensional sine- Gordon models n|, and 



monopole systems [11 



Although oscillons are not necessarily exact solutions, they live for extremely long times 
- in some cases so long that sophisticated analytic arguments are needed to decide their 
stability 0, k| because numerical experiments never see their decay. For many physics 



applications, the distinction between an infinite and a very long lifetime is irrelevant. As 
long as the object's lifetime is significantly larger than the natural time scales of the problem, 
it can have significant effects on the dynamics of the theory. 

In this Letter we display an oscillon in the gauged SU(2) Higgs model, which is the 
Standard Model of the weak interactions neglecting electromagnetism and fermions. We set 
the Higgs mass to be twice the gauge boson mass. We work in the spherical ansatz and 
numerically evolve the classical equations of motion. Even after extensive runs, exceeding 
times of 50,000 in natural units, we have never seen this oscillon decay. 



THE MODEL 



We consider an SU (2) gauge theory coupled to a doublet Higgs in 3 + 1 dimensions. The 
Lagrangian is 



--TrF^F^ + ^Tr (D^) ] - - (Tr $ f $ - v 2 )' 



where 

= d fl A u - d u A^ - ig[A^ A v \ , = (fy - igAJQ , A^ = A«a a /2 , 

and we have defined the 2x2 matrix $ to represent the Higgs doublet (p by 



ip 2 ip x 



(1) 



(2) 



(3) 



We follow the conventions of [l4|, except we use the metric ds 2 = dt 2 — dx 2 . 

The spherical ansatz [l5| is given by expressing the gauge field A^ and the Higgs field $ 
in terms of six real functions ao(r,t), ai(r,t), a(r,t), 7(r, t), fi(r, t) and v{r, t): 



A (x,t) = —a {r,t)cr ■ x 

A(x,t) = — 

&(x,t) = -[//(r, t) + iu(r, t)cr ■ x] 
9 



/ i\ . oi(r,t) 7(r,t) 
ai(r, t)a ■ xx^ H (t7j - a ■ xxj H e ijk Xja k 



(4) 
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where x is the unit three-vector in the radial direction and cr are the Pauli matrices. For 
the fields of the full theory to be regular at the origin, ao, a, a\ — a/r, j/r and v must 
vanish as r — > 0. The theory reduced to this ansatz has a residual U(l) gauge invariance 
consisting of gauge transformations of the form exp[ifl(r,t)cr ■ x/2] with Q(0, t) = 0. 
In the spherical ansatz we obtain the reduced Lagrangian density 



£(r,t) = — 



\r 2 rU» + (D^xYD.x + r 2 {D^yD,4> - ± (\ x \ 2 - l) : 



i(|xr + l)|0| 2 -Re(^0 2 )-4 r2 
^ 9 



2 2\ 2 



(5) 



where the indices now run over and 1 and 



///' d^ci u <) u (i p , 



X = a + i(7 - 1) , 
DpX = {9^ - ia^x , and = 



fJL + %V . 



(6) 



Under the reduced U(l) gauge invariance, the complex scalar fields X an d 4> have charges 
of 1 and 1/2 respectively, is the gauge field, f^ u is the field strength, and is the 
covariant derivative. The indices are raised and lowered with the 1 + 1 dimensional metric 
ds 2 = dt 2 - dr 2 . 

The equations of motion for the reduced theory are 



^ 2 + ^(lxl 2 -D + ^l 



^r 2 J D M + -(| X | 2 + l) + 



2A 
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12 9 2 v 2 



X = ~2 C 



<p = ix<p*. 



(7) 



Although it is described by six fields, there are only four independent degrees of freedom 
in the theory, consisting of the three VT-bosons, each radially polarized, and the massive 
Higgs. The remaining degrees of freedom are gauge artifacts. We may fix the gauge by 
setting ao(r, t) = everywhere in space and time, and then applying a time-independent 
gauge transformation to set ai(r, t = 0) = 0. For each of the four fields fi, u, a, and 7, we 
must specify the profile at t = as a function of r and the time derivative of the profile 
at t = as a function of r. The time derivative of ai(r,t) at t = is then determined by 
imposing Gauss's Law (the first equation in (J2J) with index v — 0). Configurations obeying 
Gauss's Law at the initial time will then obey it for all times. Thus we have fully specified 
the configuration by providing initial value data for X an d 4> only, reflecting the four real 
degrees of freedom of the model. 



NUMERICAL SETUP 

In order to accurately simulate the extremely long lifetimes of oscillons, we require numer- 
ical techniques that are highly stable. We discretize the system at the level of the Lagrangian 
in Eq. (jSJ). We choose a fixed radial lattice spacing Ar, placing the scalar fields at the sites 
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of the lattice and the gauge field a% on the links. Thus each of the scalar fields £ (r, t) is 
replaced by a set of functions defined at each lattice point. We use a first-order 

lattice gauge-covariant derivative given by the replacement 



£ {n+1} (t) exp[-iea 



{n+\} 



(f)Ar]-£<»>(f) 



D r £(r,t) 



i 



(8) 



Ar 



where e is the charge of the scalar field and n = r j At labels the lattice point corresponding to 
radius r. We then vary the spatially discretized Lagrangian to obtain second-order accurate 
lattice equations of motion. In the limit of continuum time evolution, this system has exact 
conservation of energy and exact gauge invariance, so we can monitor these quantities to 
detect numerical errors and determine whether our time steps are short enough. (They 
do not, however, tell us whether our spatial grid is fine enough, because even for a very 
coarse grid our system conserves energy and is gauge invariant.) In all our simulations, 
these invariants hold to an accuracy of roughly one part in 10 6 , and we see no numerical 
instabilities even after extremely long runs. 

Since we have a Q (r,t) = 0, the covariant time derivative coincides with the ordinary time 
derivative. Our time evolution is simply second-order differencing with a fixed time step 
At = Ar/2, in which each subsequent time step is computed from the previous two. This 
approach is stable for At < Ar. We note that to monitor the energy conservation and gauge 
invariance of the solution, it is necessary to compute first-order time derivatives. These 
derivatives are computed by subtracting the result at to from the result at to + 2 At, to get 
a second order accurate result at to + At. 



We found the oscillon by picking localized initial value data with zero time derivatives, and 
letting these configurations evolve. For some initial values, the fields dissipate completely. 
However, for others, much of the energy radiates away but an oscillating localized core 
remains. If we then start from a configuration close to this core, we can find a configuration 
for which relatively little energy radiates away. 

We work in natural units with v — 1. Since we are considering only classical dynamics, 
the theory is invariant under overall rescalings of the Lagrangian density. As a result, the 
theory is completely specified by the choice of the ratio of the Higgs mass m# = vy2X to 
the H^-boson mass mw = gv/2. For generic values of this ratio, we observe configurations 
that oscillate with extremely gradual decay. For the particular case m# = 2mw, however, 
these oscillations exhibit no observable decay whatsoever. We choose g = y/2 and A = 1 
and start from the localized configuration 



(j)(r,t = 0) = ±= (l + d l e- r2/w2 ) and x(r, t = 0) = -% ( 1 + id 2 e~ r2/w2 ) , (9) 



where d\, d 2l and w parameterize the configuration. The initial time derivatives are set 
equal to zero everywhere. 

We display results for d x = —0.1, d 2 = —3.2, and w = 12, although small changes in these 
parameters give similar results. In the left panel of Figure ^ we show the energy density 
as a function of t and r. A localized, time dependent object is clearly visible. The right 
panel shows that the initial configuration sheds about a quarter of its energy and quickly 



RESULTS 
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FIG. 1: Energy density as a function of position and time (left). Total energy and energy in the 
box of size 80 (right). 

settles into the localized oscillon. In the upper left panel of Figure |2] we show the gauge 
invariant quantity \<p(r = 0,t)\. The field is clearly oscillating about the vacuum expectation 
value at = 1. The frequency of the oscillation is measured to be fn = 0.2239. A general 
characteristic of any long-lived oscillon is that it oscillates at a frequency below the frequency 
of the modes that can carry energy away [T^f. Note that the mass of the Higgs is m# = \/2, 
so the lowest frequency of any propagating mode of the Higgs field is m H /2ir = 0.2251, 
which is just above the oscillon frequency. In the upper right panel we show the gauge 
invariant quantity /oi( r — 10, t), which is oscillating with a frequency of fw = 0.1120. The 
gauge boson mass is m w = y/2/2, so the lowest frequency of these propagating modes is 
mw /27r = 0.1125, which is also just above the oscillation frequency. The field x{ r i t) has the 
same mass as foi (r, t) and oscillates at the same frequency. 

Oscillons are localized in space. Going far from the origin, where we can linearize the 
theory, for each field we expect to find a relationship of the form 2rrf = \Jm 2 — e 2 , where the 
fields approach their vacuum values like exp(— er) for large r. Unfortunately, the precision 
with which we can estimate the e for each field from our numerical data is very limited, so we 
have only checked that it agrees with the predicted value to within an order of magnitude. 

In the lower panels of Figure El we plot the two fields over a much longer time, so the 
fundamental oscillations are not visible. On this scale, however, we see beats, at the same 
frequency in all fields. For this configuration the beat frequency is about 5 x 10~ 4 , but it 
can vary substantially in response to variations of the initial parameters. 

CONCLUSIONS 

In the model we have studied, we have numerically discovered an oscillon that shows 
no sign of decay after more than 14,000 oscillations. This result suggests that this object 
is extremely long-lived. If such a similarly long-lived object were to exist in the Standard 
Model, it would have important phenomenological consequences. It could play a significant 
role in the dynamics of baryogenesis by providing a mechanism for creating the necessary out- 
of-equilibrium conditions in the early universe. If this object were extraordinarily long-lived, 
it could be a dark matter candidate. (A substantial body of work has already investigated 
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FIG. 2: \4>(r = 0,t)\ (left) and /oi(r = 10, i) (right) as functions of time for d\ = —0.1 and 
d 2 = -3.2. 

such questions for Q-balls fl7l fisL fToj| . ) 

However, a number of questions must be resolved in order to understand the phenomeno- 
logical implications, if any. Although one might expect objects with spherical symmetry to 
be energetically favored, it is possible that our oscillon could decay by coupling to nonspher- 
ical deformations, which we would not see because of our restriction to the spherical ansatz. 
Restoring the photon coupling would add another possible decay mechanism. The oscillon 
could also decay by coupling to fermions, and such decays might be of interest in the early 
universe. It would be interesting to see if quantum effects modify these conclusions, as they 
do for small Q-balls |2(|. Finally, and perhaps most importantly, we have only observed 
oscillons when m H = 2m w , so it is possible that their unnaturally long lifetime is a conse- 
quence of this fine-tuning of parameters in our reduced theory. We imagine that in the full 
theory, the condition needed for the existence of the oscillon would be slightly modified by 
the splitting between the W and Z masses. We note that a Higgs mass near 2m\y is in the 
middle of the discovery window of the LHC. 
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